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ABSTRACT 


The angular size of the broad line region (BLR) of the nearby active galactic nucleus (AGN) NGC 3783 has been spatially resolved by recent 


= observations with VLTI/GRAVITY. A reverberation mapping (RM) campaign has also recently obtained high quality light curves and measured the 
linear size of the BLR in a way that is complementary to the GRAVITY measurement. The size and kinematics of the BLR can be better constrained 
N by a joint analysis that combines both GRAVITY and RM data. This, in turn, allows us to obtain the mass of the supermassive black hole in 
NGC3783 with an accuracy that is about a factor of two better than that inferred from GRAVITY data alone. We derive Mgy = 2.54020 x 107 Mo. 


r—4Finally, and perhaps most notably, we are able to measure a geometric distance to NGC 3783 of 39.9 


+14.5 


^1» Mpc. We are able to test the robustness 


of the BLR-based geometric distance with measurements based on the Tully-Fisher relation and other indirect methods. We find the geometric 
distance is consistent with other methods within their scatter. We explore the potential of BLR-based geometric distances to directly constrain the 
Hubble constant, Ho, and identify differential phase uncertainties as the current dominant limitation to the Ho measurement precision for individual 


Ea sources. 
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1. Introduction 


astro 


—— Trigonometry is the basis of distance measurements. The paral- 
lax method uses the motion of the Earth around the Sun to mea- 
sure the angular displacement of a nearby star (Bessel 1838). 
From the Hipparcos satellite of the European Space Agency 

\© (ESA; ESA 1997) to the recent Gaia mission (Gaia Collabora- 

CN tion et al. 2016), the parallaxes of more than 1 billion stars in 

SF the Milky Way have now been measured (Gaia Collaboration 

73 et al. 2018). The geometric method can also be applied to any 

* object, as long as both its physical size (R) and angular size 
© (©) are measurable, as D = R/O. For a distant extragalactic 
v4 target, the measured geometric distance is the angular diameter 

C distance (D4), including cosmological expansion. In contrast to 

—* "standard candles,” such as pulsating stars (e.g., Bhardwaj 2020) 

«= and Type Ia supernovae (SNe Ia; e.g., Phillips 1993; Riess et al. 

1996), which measure the luminosity distance (Dz), the geomet- 

= ric method does not rely on the calibration based on the so-called 
C3 distance ladder (e.g., Riess et al. 2009, 2021). 

However, objects to which the geometric method can be ap- 

plied are usually rare as it is difficult to measure both R and © 


* GRAVITY is developed in a collaboration by the Max Planck 
Institute for extraterrestrial Physics, LESIA of Observatoire de 
Paris/Université PSL/CNRS/Sorbonne Université/Université de Paris 
and IPAG of Université Grenoble Alpes / CNRS, the Max Planck Insti- 
tute for Astronomy, the University of Cologne, the CENTRA - Centro 
de Astrofisica e Gravitacáo, and the European Southern Observatory. 
** Corresponding author: J. Shangguan 
e-mail: shangguan ? mpe.mpg.de 


for the same target. The distance to the supermassive black hole 
(BH) at the center of the Milky Way is measured to a « 0.596 
uncertainty level, based on the 27-year astrometric and spectro- 
scopic monitoring of the 16-year orbital motion of the star S2 
(Gravity Collaboration et al. 2019). For detached eclipsing bi- 
nary stars, the linear size of each component can be measured 
from photometric and spectroscopic monitoring, while the an- 
gular size of each star can be derived from its empirical relation 
with the color of the star (Lacy 1977). This method has been 
used to measure the distance to the Large Magellanic Cloud with 
high accuracy (Pietrzyfiski et al. 2019). By monitoring the Kep- 
lerian motion of water maser-emitting gas with very long base- 
line interferometry (VLBI) observations, one can measure the 
geometric distance to nearby Type-2 Seyfert galaxies that host 
an observable megamaser disk (Herrnstein et al. 1999; Braatz 
et al. 2010). Baryonic acoustic oscillations can also provide D4 
using the clustering of galaxies at a certain redshift range (e.g., 
Eisenstein et al. 2005; Anderson et al. 2014). The broad line re- 
gion (BLR) of active galactic nuclei (AGNs) has also been pro- 
posed as a probe of geometric distance (Elvis & Karovska 2002). 
The linear size of the BLR can be measured by the reverberation 
mapping (RM) technique (Peterson 2014), while the angular size 
of the BLR can be measured by near-infrared (NIR) interferom- 
etry (Petrov et al. 2001; Woillez et al. 2004). Unlike detached 
eclipsing binaries and megamaser systems, which are difficult to 
discover, AGNs are luminous sources that are commonly found 
locally as well as out to high redshifts, even at z > 6. A simi- 
lar approach has been applied to NGC 4151 by resolving the hot 
dust continuum emission instead of the BLR (Hónig et al. 2014). 
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Recently, GRAVITY, the second generation Very Large Tele- 
scope Interferometer (VLTT) instrument, spatially resolved the 
BLR of 3C 273 for the first time using the spectroastrometry 
(SA) technique (Bailey 1998). GRAVITY combines all four of 
the 8 m Unit Telescope (UT) beams to yield six simultaneous 
baselines (Gravity Collaboration et al. 2017). Gravity Collab- 
oration et al. (2018) reported the mean radius of the BLR as 
46 + 10 uas and a BH mass of (2.6 + 1.1) x 108 Mo, which is 
fully consistent with that measured by Zhang et al. (2019) us- 
ing 10 year RM data. Wang et al. (2020) conducted the first 
joint analysis (hereafter, SARM, as the combination of SA and 
RM) and derived an angular diameter distance for 3C 273 of 


Da = 551.5213 Mpc, corresponding to a Hubble constant 


Ho = 71.5+1}? kms! Mpc™!. 

The Hubble constant has been measured to an accuracy of a 
few percent in the low-z Universe with various methods. Using 
the distance ladder, the SHOES (supernovae HO for the equation 
of state) project recently derived Hy = 73.5 + 1.4km s7! Mpc”! 
(Reid et al. 2019). This value is consistent with other measure- 
ments independent of the distance ladder, such as megamaser 
observations (e.g., Pesce et al. 2020) and observations of grav- 
itationally lensed quasars (e.g., Wong et al. 2020) based on the 
so-called time-delay distance. In comparison, Ho inferred from 
conditions in the early Universe tends to be substantially smaller 
than the low-z values (4—6 c significance level, Riess 2020; how- 
ever, see Freedman et al. 2019). For instance, Planck Collabo- 
ration et al. (2020b) derive 67.4 + 0.5 kms^! Mpc™! based on 
cosmic microwave background observations. The high luminos- 
ity and remarkably uniform spectral properties of AGNs have 
motivated many attempts to use them as standard candles (Bald- 
win 1977; Collier et al. 1999; Elvis & Karovska 2002; Watson 
et al. 2011; Wang et al. 2013; Honig 2014; La Franca et al. 2014; 
Risaliti & Lusso 2015). In particular, the good correlation be- 
tween the BLR radius and continuum luminosity (e.g., Kaspi 
et al. 2000; Bentz et al. 2013) makes it promising to probe the 
Dy, based on the RM measurements (e.g., Watson et al. 2011; 
Wang et al. 2013). However, recent studies reveal significant de- 
viations from the R-L relation, primarily driven by the increase 
in the Eddington ratio (Du et al. 2016, 2018; Du & Wang 2019; 
Martínez-Aldama et al. 2019; Dalla Bontà et al. 2020; however, 
see Grier et al. 2017; Fonseca Alvarez et al. 2020). It is therefore 
of great importance to explore the power of the SARM method, 
namely to study the BLR structure in detail and to measure the 
BLR-based geometric distance with the goal to independently 
test the Ho tension in the future. 

With z ~ 0.009730 (Theureau et al. 1998) and a much shorter 
time lag than that of 3C 273, NGC 3783 provides the best oppor- 
tunity to compare the geometric distance derived with SARM 
to other independent measurements. A new RM campaign of 
NGC 3783 was reported by Bentz et al. (2021) with a time lag of 
about 10 days, consistent with previous measurements (Onken 
& Peterson 2002). Gravity Collaboration et al. (2021, hereafter, 
GC21) reported a BLR mean angular radius of about 70 uas, and 
the RM-measured time lag can be reproduced with the measured 
continuum light curve and the best-fit BLR model inferred only 
from GRAVITY data at an assumed D4 = 38.5 Mpc. In this pa- 
per, we construct a Bayesian model to fit the GRAVITY and RM 
data simultaneously (Section 3). The inferred BLR model is en- 
tirely consistent with the results of GC21 (Section 4). We find 
that the inferred D4 of NGC 3783 is fully consistent with dis- 
tances measured using the Tully-Fisher relation and other indi- 
rect methods (Section 5). The application of the SARM method 
to Ho is discussed in Section 5.4, together with improvements 
that may come with the ongoing upgrade of the GRAVITY in- 
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strument. This work adopts the following fiducial parameters for 
a ACDM cosmology: Qm = 0.3, Qa = 0.7, and Ho = 70 km s7! 
Mpc^!, unless otherwise specified. 


2. GRAVITY interferometer and reverberation 
mapping observations 


GRAVITY interferometric data were collected over 3 years from 
2018 to 2020 through a series of Open Time programs and an 
ESO Large Programme! with the aim to measure the size of the 
BLR and the mass of the central BH. Details of the data reduc- 
tion and analysis can be found in GC21. We only briefly sum- 
marize the main points here. Phase referenced with its hot dust 
continuum emission, NGC 3783 was observed with MEDIUM 
spectral resolution (R ~ 500) in the science channel and com- 
bined polarization. After the pipeline reduction, we exclude data 
with poorer performance in terms of phase reference, keeping 
those with fringe tracking ratio > 80%. We find that the Bry 
profiles in the GRAVITY spectra are consistent considering the 
calibration uncertainty, which indicates that the BLR did not 
change significantly between observations. We therefore aver- 
age the differential phase of all of the data into three epochs ac- 
cording to their uv coordinates. The differential phase signal of 
the BLR is dominated by the so-called continuum phase, indi- 
cating that the center of the BLR is offset from the photocenter 
of the hot dust continuum emission. The reconstructed image of 
the NGC 3783 continuum emission displays a secondary com- 
ponent, which is the main cause of the continuum phase signal, 
although the asymmetry of the primary continuum emission can 
also contribute to the continuum phase (Gravity Collaboration 
et al. 2020, hereafter, GC20). Following GC21, for our analysis 
we adopt the differential phase data after subtracting the contin- 
uum phase that is calculated from the reconstructed continuum 
image. There could be still some residual continuum phase due 
to the uncertainty of the spatial origin of the image reconstruc- 
tion, so we still allow the BLR to be offset from the continuum 
reference center in the modeling (Section 3.1). The amplitude of 
the differential phase due to BLR rotation is about 0.2? (Figure 
8(b) of GC21). 

Following GC21, for the model inference we adopt the Bry 
profile from our R « 4000 adaptive optics observations obtained 
on 20 April 2019 using SINFONI (Eisenhauer et al. 2003; Bon- 
net et al. 2004). The nuclear spectrum is extracted using a cir- 
cular aperture centered on the peak of the continuum and a di- 
ameter of about 0.1 arcsec, which is larger than the field of view 
of GRAVITY ~ 60 milliarcsec. This results in higher narrow 
Bry emission in the SINFONI spectrum; meanwhile, the broad 
Bry profile from the SINFONI spectrum is consistent with those 
from GRAVITY observations. The high spectral resolution of 
the former guarantees a robust decomposition of the narrow-line 
component, while the latter suffer from low spectral resolution 
and systematic uncertainties due to the calibration. Therefore, 
we prefer to use the SINFONI data in the BLR modeling. 

The RM data were collected in the first half of 2020. The de- 
tails of the observations and analyses are reported in Bentz et al. 
(2021). Briefly, the photometric and spectroscopic monitoring 
was conducted with the Las Cumbres Observatory global tele- 
scope (LCOGT) network. A total of 209 V-band images and 50 
spectra were obtained. The [O 111] 444959,5007 doublet region 
was used to calibrate the relative flux of the spectra, resulting in 


! Observations were made using the ESO Telescopes at the La Silla 
Paranal Observatory, program IDs 0100.B-0582, 0101.B-0255, 0102.B- 
0667, 2102.B-5053, and 1103.B-0626. 
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a ~ 2% accuracy for relative spectrophotometry of H£ through- 
out the monitoring. The emission line light curves were derived 
from the calibrated spectra by integrating the emission lines with 
the local continuum fitted and subtracted. The continuum light 
curve at rest-frame 5100 À is also measured from the spectra 
and it is used to cross-calibrate the V-band light curve. Both of 
them are merged into the final continuum light curve. The RMS 
spectrum shows that Hf, He 11 44686, Hy, and Hô are variable. 
A time lag of about 10 days is derived from the well calibrated 
light curves of the continuum and Hf line. 

Previous studies have found that the higher order Balmer 
lines tend to show shorter time lags (e.g., Kaspi et al. 2000; Bentz 
et al. 2010). High-ionization lines such as He 11 44686 also show 
smaller lags than that of low-ionization lines such as Balmer 
lines (e.g., Clavel et al. 1991; Peterson & Wandel 1999; Bentz 
et al. 2010; Grier et al. 2013; Fausnaugh et al. 2017; Williams 
et al. 2020). Photoionization models (Netzer 1975; Rees et al. 
1989; Baldwin et al. 1995; Korista & Goad 2004) provide phys- 
ical explanations for the radial stratification of the BLR. The 
higher order Balmer lines, with lower optical depth, are expected 
to be more efficiently emitted by gas with higher density and 
closer to the central BH, resulting in a shorter time lag and higher 
"responsivity" (see Korista & Goad 2004). However, the optical 
depth of H8 is z 10? based on typical BLR models (e.g., Net- 
zer 2020). Itis very challenging for photoionization models such 
as CLOUDY (Ferland et al. 2017) to theoretically calculate the 
Hydrogen line emission. 

Our joint analysis assumes that RM and GRAVITY probe 
the same regions of the BLR, which is encouraged by the con- 
sistent size measured from the H6 time lag and GRAVITY mea- 
surements of the Bry line. Wang et al. (2020) conducted a joint 
analysis using the H£ light curve and GRAVITY measurements 
of the Pao line. The choice is mainly limited by the available 
measurements, although Zhang et al. (2019) show that the time 
lags of H8 and Hy for 3C 273 are consistent. It is still unclear 
however how well the BLR sizes of Hydrogen lines in the op- 
tical and NIR are consistent with each other. Unfortunately, the 
Hy and Hé light curves of NGC 3783 are not robustly calibrated, 
because they are far from the reference [O i11] lines. Their time 
lags are 2-4 times smaller than that of HG, while their FWHMs 
are also smaller than that of Hf. This strongly suggests that the 
lags of Hy and Hó are not physically robust. Bentz et al. (2021) 
also report a ~ 2-o detection of the lag of Herr about 2 days; 
however, the discrepancy of BLR sizes between H8 and Her 
lines are not unexpected. We therefore only adopt the H£ light 
curve together with the continuum light curve in our joint anal- 
ysis. One of our main goals in this work is to test whether HS 
and Bry BLR radii are consistent by comparing our geometric 
distance with other independent distance measurements. 


3. SARM joint analysis 
3.1. BLR model and spectroastrometry 


Our BLR model has been introduced in detail by Gravity Col- 
laboration et al. (2018) and GC20. Here, we only provide a 
brief description of the model that is necessary for this work. 
The model was first developed by Pancoast et al. (2014a) with 
the original purpose to model velocity resolved RM data. The 
BLR is assumed to consist of a large number of non-interacting 
clouds, whose motion is governed only by the gravity of a cen- 
tral BH with mass Mpy. A shifted gamma distribution, r = 
Rs + F Rgig * gd — Fg Rgxr, is used to describe the radial 
distribution of the clouds, where Rs is the Schwarzschild radius, 


g = p(x|1/8?, 1) is drawn randomly from a Gamma distribution, 
p(ala, b) = x*-le7*/^ (T(a) b^), and I (a) is the gamma function. 
The shape parameter, B, controls the radial profile to be Gaussian 
(0 < B < 1), exponential (8 = 1), or heavy-tailed (1 < B < 2). 
The weighted mean cloud radius, Rp; p, and fractional inner ra- 
dius, F = Rin/Rpır, are fitted. Clouds are then randomly dis- 
tributed in a disk with the angular thickness 6), which ranges 
from 0? (thin disk) to 90? (sphere). y controls the concentra- 
tion of the cloud distribution toward the edge of the disk (Equa- 
tion (4) of GC20). The weight of the emission from the clouds is 
controlled by «. With respect to the observer, the near side clouds 
have higher weight if x > 0, while the far side clouds have higher 
weight otherwise (Equation (5) of GC20). The fractional differ- 
ence in the number of clouds above and below the mid-plane is 
controlled by £. Setting € = 0 means there are equal numbers 
of clouds each side of the mid-plane, while there are no clouds 
below the mid-plane if £ = 1. 

We define the velocity of the clouds according to their posi- 
tion and the BH mass. We draw cloud velocities from the param- 
eter space distribution of radial and tangential velocities centered 
around either the circular orbits for bound clouds or around the 
escape velocity for inflowing or outflowing clouds (Equation (6) 
of GC20). The fraction of clouds in bound elliptical orbits is 
controlled by fenip. Clouds in radial orbits (inflowing or outflow- 
ing) are allowed to be mostly bound, mainly controlled by 6. 
A random distribution is drawn with Gaussian dispersion along 
and perpendicular to the ellipse connecting circular orbit veloc- 
ity and radial escape velocity. Whether the radial motion of the 
cloud ensemble is inflowing or outflowing is controlled by a sin- 
gle parameter ffow, < 0.5 for inflowing and > 0.5 for outflowing. 
The model further considers a line-of-sight velocity dispersion to 
model the macroturbulence. As in GC21, we find the dispersion 
parameters of the Gaussian distributions in the phase space are 
not crucial to fit the data. We fix them to zero, so that the BLR 
model is the same as that adopted in GC21. 

The BLR model is rotated with inclination angle i, ranging 
from 0° (face-on) to 90° (edge-on), and position angle PA. Line- 
of-sight velocities of the clouds account for the full relativis- 
tic Doppler effect and gravitational redshift. The flux of each 
spectral channel is calculated by summing the weights of clouds 
in each velocity bin. The model line profile is scaled according 
to the maximum, fpeak. The photocenter of each channel is the 
weighted average position of all clouds in each bin. The differ- 
ential phase at wavelength J is calculated as 

fa 
l+ fa 
where f; is the line flux at wavelength 4 to a continuum level of 
unity, u is the uv coordinate of the baseline, and x, is the offset of 
BLR center from the reference center of the continuum emission, 
which can introduce the continuum phase (GC20). 


Ag, = -27 U : (XBrga — Xc), (D 


3.2. Light curve modeling 


In order to generate the model Hf light curve, we need to cal- 
culate the continuum flux reverberated by clouds with different 
time lag at a given observed time. We need, therefore, to inter- 
polate and extrapolate the continuum light curve taking into ac- 
count the measurement uncertainty. The variability of an AGN 
can be described by a damped random walk model (Kelly et al. 
2009) for which the covariance function between any two times 
tı and f» is 


t-t 
Sti, h) =o} ap[- 525. Q) 
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where c4 is the long-term standard deviation and Ty is the typi- 
cal correlated timescale of the continuum light curve. We use a 
Gaussian process to model the continuum light curve (e.g., Pan- 
coast et al. 2011, 2014a; Li et al. 2018). In the fitting, we adopt 
T4 and Ôa = og/ ra in order to relax the correlation between 
the two parameters. 

We can calculate the model H£ light curve by convolving the 
model continuum light curve l(t) (Appendix A) with a so-called 
transfer function Y(T), 


lug(t) = A [ vo L(t — 7) dr, (3) 
where A is a scaling factor. The transfer function is the normal- 
ized distribution of the time lag taking into account the weight 
of clouds of the BLR. 


3.3. Bayesian inference 


Following GC20, we fit the observed Bry profile and differential 
phase with the likelihood function of GRAVITY spectroastrom- 


etry, 
«I G-A 
oe fi 2c E M 


(gi - ài 
x | | ——— exp : (4) 
a V2na; $.i NH EE 
where f; and fi are the observed and model fluxes of the Bry 
profile; $; and ¢; are the observed and model differential phases; 
Cy; and og; are the measurement uncertainties of f; and dj re- 
spectively; and i denotes the ith channel. The measured Hf light 
curve is fitted with the likelihood function of RM, 
_ (ingi = lug; | 


LRM = 
RM 20, 


— E | (5) 
ist V27O Hpi 


where /ug; and ogg; are the ith measurements of Hf flux and 
uncertainty. Therefore, the joint likelihood function is the multi- 
plication of Equation (4) and Equation (5), Lsarm = Lsa XLRM. 

The physical parameters of the BLR model are the primary 
parameters to be inferred from the joint analysis. Parameters of 
the continuum light curve model (xs, xg, Ta, and Gg) are also 
involved in the fitting. x, and x, are the deviation of the con- 
tinuum light curve fluxes and their long-term average from the 
maximum a posteriori conditioned by the observed light curve 
using the Gaussian process regression (Appendix A). Moreover, 
x, consists of 200 parameters because it is important to densely 
sample the model continuum light curve. It is however worth 
emphasizing that x, is well constrained by the prior information 
of the observed continuum data and it only enters the likelihood 
function via Equation (3). We tested the fitting using 100—300 
points for x,. While there is no significant difference in the fit- 
ting results, we find that the reconstructed continuum light curve 
with 100 points does not capture some features of the observed 
data in some densely sampled region. The goodness of the fit- 
ting does not improve when we adopt 300 points. Therefore, for 
the sake of computation power, we adopted 200 points in our 
analysis. With 221 free parameters in total, it is very challenging 
to sample the parameter space. We utilized the diffusive nested 
sampling code CDNest (Li 2018) to do so. Diffusive nested sam- 
pling (Brewer et al. 2011) has been shown to be effective for 
fiting RM data with a BLR model (e.g., Pancoast et al. 2014a,b; 
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Parameters Model Inference 
Oni n (uas) qim 
Rere (Id) 16.2775 
D4 (Mpc) 390755. 
log (Mau /Mo) 74098 


Table 1. Inferred median of posterior sample and central 68% credible 
interval for the key parameters of the analysis. Opie is the angular ra- 
dius of the BLR, which is mainly constrained by the differential phase 
of GRAVITY data. Rgig is the physical radius of the BLR, which is 
primarily constrained by RM data. D; is the angular diameter distance. 
Mbegu is the BH mass. 


Li et al. 2018), as well as for the joint analysis of 3C 273 (Wang 
et al. 2020). 

We find that Gq and tg are easily biased in the joint fitting 
if no informative prior is used. The covariance model (Equa- 
tion (2)) should be able to describe the continuum light curve 
well independently. Therefore, we opted to optimize the likeli- 
hood function of the continuum light curve data given the covari- 
ance model to constrain Ôq and r4 in advance of the joint fitting, 


Be 1 (- (le - Eq)’ Cur (I - 6 
O ¥Qny"Cul 2 
where Cu = S(t1,t1) + e'I is the model covariance ma- 


trix of the measured continuum light curve (le) and uncer- 
tainties (Fe); and the scalar q is the long-term average of the 
light curve (see Appendix A for more details). The flat pri- 
ors with (—2, 2) and (—1, 3) for log (ta/day) and log (64/ day) 
are wide enough for this purpose. We are able to obtain good 


constraints of log(ra/day) = 1.75:087 and log (@a/+day) = 


-0.37 
—0.75:0:2. Therefore, we adopted the priors Norm(1.75, 0.5?) 


and Norm(—0.75, 0.037) for Ta and ĉa, respectively in logarith- 
mic scale in the joint fitting. Unlike 3C 273 (Wang et al. 2020; Li 
et al. 2020), we do not find a long-term trend in the light curve of 
NGC 3783. Therefore, a simple constant g is enough to describe 
the long-term average of the light curve. We emphasize that the 
parameters that we are mainly interested in are the parameters of 
the BLR model as well as D4, while the parameters of the light 
curves are considered as nuisance parameters. 


4. BLR model inference 


The inferred model parameters from the joint analysis are en- 
tirely consistent with what we have obtained from fitting GRAV- 
ITY data only (GC21). Table 1 presents the key parameters of 
this work, while the full model parameters are reported in Ta- 
ble B.1. All of the parameters are consistent within their 2-c 
uncertainty levels. One difference in our approach here is that, 
instead of calculating the maximum a posteriori as GC21, we 
simply report the median of the marginalized posterior distribu- 
tion, because the dimension of the parameter space (> 200) is 
prohibitively high to robustly calculate the maximum a posteri- 
ori. Nevertheless, we confirm that the median and maximum a 
posteriori of the posterior sampling using only GRAVITY data 
always show difference within 1-c. The fitting results are dis- 
cussed in detail in Appendix B. We would like to highlight that 
the uncertainties of several key parameters, in particular Mpy, 
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Fig. 1. Best-fit line profile and differential phase signal of the BLR are 
displayed with red curves. The Bry profile from SINFONI spectrum is 
in black. The averaged differential phase data of UT4-UT2, UT4-UT1, 
and UT3-UT1 (see Figure B.1) with the continuum phase subtracted are 
displayed in blue. 


are reduced by a factor of ~ 2 after adding the RM light curves 
in the fitting, which indicates that RM data help to constrain 
the model in a consistent manner with GRAVITY data (see also 
Wang et al. 2020). Our inferred BH mass, 2.54*099 x 107 Mo, is 
fully consistent with the value based on RM-only data by Bentz 
et al. (2021). Our uncertainty is about a factor of two larger than 
their statistical uncertainty; but one should bear in mind that this 
does not include the ~ 0.3 dex uncertainty due to the virial fac- 
tor, which is the primary uncertainty of integrated RM (e.g., Ho 
& Kim 2014; Batiste et al. 2017). Our best-fit model favors mod- 
erate inflow. This is expected because both the GRAVITY data 
and RM data indicate a preference for gas inflow individually. 
Bentz et al. (2021) derived the time lags of H£ in 5 velocity bins. 
The time lag profile across the line is slightly asymmetric with 
the longest wavelength showing the shortest time lag, indicating 
that the gas motion of the BLR is a combination of rotation and 
inflow (however, see Mangham et al. 2019). 


The line profile and differential phase signal based on the in- 
ferred median parameters are plotted in Figure 1. In this figure, 
the differential phase is averaged for UT4-UT2, UT4—UTI, and 
UT3-UT1. The continuum phase (Section 2) is subtracted based 
on the median offset from the data of each baseline before the 
average. The data averaged phase displays moderate asymmetry, 
the excess in data points with respect to the best-fit model seen 
on the blue side of the “S” shape, which can be explained by the 
uncertainty of those offsets. The reconstructed continuum light 
curves and the reverberated Hf light curves are shown with the 
observed data in Figure 2. The detailed features of the contin- 
uum light curve are captured by the model, while the H£ light 
curve is reasonably well fitted. We infer a correlation timescale 
Tq & 87 day, slightly longer than the center of the prior distribu- 
tion. The mean time lag of the BLR model, based on the median 


inferred parameters, is 14.173 day, which is close to the light- 


weighted radius of the BLR, 16.2*> Id (Figure 3). However, the 


-3.5 
centroid of the cross correlation function (CCF) of the model 
continuum and Hf light curves is Tent = 8.2 day. Following 


Bentz et al. (2021), we derive Tent as the first moment of the CCF 
above 0.8 of the peak of the CCF (Koratkar & Gaskell 1991). 
Our Tee is slightly lower than, but within 2-c to, that reported 
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HJD-2450000 (days) 


Fig. 2. Fitting results of (a) continuum and (b) Hf light curves. We 
randomly selected 50 reconstructed continuum light curves and the re- 
verberated Hg light curves from the posterior sample and display them 
in gray. The median light curves are in red. The data points are plotted 
in blue. 
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Fig. 3. Posterior probability distribution of the BLR radius, BH mass, 
and angular diameter distance of the joint analysis. The dashed lines 
indicate the 16%, 50%, and 84% percentiles of the of the posterior dis- 
tributions. The contours indicate lo, 1.50, and 20. 


by Bentz et al. (2021). This is consistent with the conclusion in 
GC21: The measured time lag from the CCF underestimates the 
BLR radius of NGC 3783 because the BLR radial distribution is 
strongly heavy-tailed (large 8; Figure 4). Nevertheless, detailed 
comparisons of BLR structures of Bry and H£ lines, including 
the velocity-resolved RM, are needed in the future. We discuss 
the fitting results considering an outer truncation of the BLR as 
well as the effect of nonlinear response of the continuum light 
curve in Appendix B. The primary model introduced above is 


Article number, page 5 of 14 


A&A proofs: manuscript no. ms 


VLos (103 km s71) 
-3 0 3 


A Dec. (uas) 


A R. A. (uas) 


Tmean = 14.1 days 
Tcent = 8.2 days 


Y( 


or- 
=] 
[e] 


0 50 60 70 
Time lag (day) 


Fig. 4. Transfer function (Equation (3)) of the best-fit model peaks at the 
time lag close to Teen, measured from the CCF method. But the mean 
lag (Tmean) is longer because the transfer function has a long tail, which 
extends beyond the limit of the plot. The inset on the right displays 
the cloud distribution of the best-fit BLR model. The color code is the 
light-of-sight velocity. 


preferred to interpret the data. While no disagreement is found 
statistically significant by introducing more physical constraints 
to the model, we caution that the inferred distance may be biased 
by the BLR model assumptions. 


In order to assess the reliability of the uncertainties on the 
derived parameters, we considered the fitting process itself as 
well as exploring how sensitive they are to the light curve. It is 
possible to adjust the temperature (T > 1), by which the log- 
arithmic likelihood is divided, with CDNest in order to enlarge 
the uncertainty of the data when the model is not flexible enough 
to fit the data (Li 2018; Brewer et al. 2011). We are always able 
to find the clear peak of the posterior weights as a function of 
the prior volume with T = 1, indicating that the fitting has prop- 
erly converged. The fitting results with T > 1 are consistent with 
those with T = 1, although the uncertainties increase. Therefore, 
we always report the fitting results obtained using 7 = 1. The 
measured uncertainty of the Hf light curve is typically ~ 1% of 
the flux, which is lower than the typical uncertainty based on the 
intercalibration using the [O ri1] doublet (Bentz et al. 2021). We 
fit the data with the Hf light curve uncertainties increased to be 
2% of the flux if they are smaller than the latter, and find that the 
results are fully consistent with those using the measured Hf un- 
certainties. The uncertainty of Mpy reaches 0.2 dex, while that 
of D4 barely increases as it is dominated by the phase data. In 
order to further test whether the joint analysis is sensitive to the 
light curve measurement, we measured the Hf light curve by de- 
composing the spectra with a relatively wide wavelength range 
(see Appendix B.3 for more details). Using the decomposition 
method, we measured the light curve with an approximatively 
2096 lower variation amplitude. We are able to obtain consis- 
tent fitting results once the nonlinear responsivity is taken into 
account in the joint analysis. 
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5. Distance and peculiar velocity of NGC 3783 


Our joint analysis infers that the angular diameter distance of 


NGC 3783 is 39.94143 Mpc. NGC 3783 has an observed red- 


shift of 0.009730 and heliocentric velocity v; ~ 2917 km s7! 
(Theureau et al. 1998). In this section, we derive the peculiar ve- 
locity of NGC 3783, compare our measured distance with other 
direct and indirect methods, and discuss the uncertainty of our 
measurements in detail. 


5.1. Peculiar velocity of NGC 3783 


To correct for the motion of the Sun in the Milky Way and the 
peculiar motion of the Milky Way, we calculate the velocity of 
NGC 3783 in the frame of the Local Sheet (Tully et al. 2008; 
Kourkchi et al. 2020), 


Vis = Vp — 26 cos Lcos b + 317 sin Lcos b — 8 sin b (7) 
= 2628kms"!, 


where (/,b) are the Galactic longitude and latitude. The Local 
Sheet reference frame is a variant of the Local Group rest frame 
(Yahil et al. 1977; Karachentsev & Makarov 1996). Tully et al. 
(2008) advocated that the Local Sheet is preferable to the Local 
Group because it is more stable. For comparison, NGC 3783 has 
a Local Group velocity 2627 km s™!, according to NASA/IPAC 
Extragalactic Database (NED) velocity calculator. The redshift 
corrected for the peculiar motion of the observing frame is 
Zi = Vis/¢ & 0.008767, which may still deviate from the cos- 
mological redshift in the Hubble flow due to the peculiar motion 
of NGC 3783. Assuming (Qm, O4, Ho) = (0.3,0.7, 70), we can 
infer a cosmological redshift z — 0.0094*0003: according to the 
measured D4. The peculiar velocity can be calculated (Davis & 
Scrimgeour 2014), 


Zis — Z 843 — 
Wee eq. = -195* T5, kms" !. (8) 


This is very close to, albeit with large uncertainties, the pecu- 
liar velocity estimated with the Cosmicflows-3 distance-velocity 
calculator (see below), —182 km s^!. We also estimate the pecu- 
liar velocity using the 6dF galaxy redshift survey peculiar veloc- 
ity map (Springob et al. 2014). We find 11 galaxies in the 6dF 
peculiar velocity catalog within 847! Mpc centered on the po- 
sition of NGC 3783. We averaged their peculiar velocities and 
uncertainties with equal weight. The estimated peculiar veloc- 
ity of NGC 3783, —158 +43 km s^!, is consistent with the other 
estimates above.? 


5.2. Comparison with other distance measures 


Besides our measured D4, the distance of NGC 3783 can be 
measured "directly" with the Tully-Fisher relation (Tully & 
Fisher 1977) on the one hand and estimated "indirectly" with 
various methods relying on the large-scale structure and velocity 


2 We adopt the recession velocity, 3197 km s™! in the CMB frame, of 
the galaxy group that comprises NGC 3783 (Kourkchi & Tully 2017) 
in order to obtain the 3-D supergalactic coordinate of NGC 3783. We 
adopt the typical smoothing scale, 847! Mpc, to find the galaxies sharing 
the same large-scale structure with NGC 3783. We find it makes little 
difference if we weight the peculiar velocity according to their separa- 
tion or not, so we average them with equal weight for simplicity. The 
peculiar velocity has been converted from the CMB frame (Planck Col- 
laboration et al. 2020a) to the Local Sheet frame in order to be compared 
with other estimates. 
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Distance (Mpc) Method description 
39.9+143 SARM (this work) 
49.8 + 19.6 Tully-Fisher relation (Robinson et al. 2021) 
38.5 + 14.2 Tully-Fisher relation (Tully & Fisher 1988) 
35.1 Cosmicflows-3 (EDD; Kourkchi et al. 2020) 
37.9 NED (Virgo + GA + Shapley) 
42.1 Galaxy group (Kourkchi & Tully 2017) 


Table 2. Distance of NGC 3783 measured with different methods. The 
SARM and Tully-Fisher relation are the direct methods, while the other 
methods are indirect. The distance uncertainties of the indirect methods 
are discussed in the text whenever they can be estimated. 


field of the local universe on the other. Strictly speaking, these 
methods, as discussed below, provide the luminosity distance. 
Since the redshift of NGC 3783 is too low for there to be a sig- 
nificant difference between angular diameter distance and lumi- 
nosity distance (< 2%), we do not distinguish the two in the fol- 
lowing discussion. The direct distance measurements and indi- 
rect estimates that we discuss below are summarized in Table 2. 

The distance of NGC 3783 has been measured as 38.5 + 
14.2 Mpc by Tully & Fisher (1988) using the Tully-Fisher re- 
lation. This measurement was based on the magnitude of the 
galaxy without removing the nuclear emission from the BH ac- 
cretion. This will bias the distance toward a lower value. Robin- 
son et al. (2021) recently derived the distance of NGC 3783 to be 
49.8 + 19.6 Mpc. They derived the maximum rotational velocity 
measured from the newly observed H 1 spectrum and measured 
the multiband optical/NIR magnitudes of the galaxy after care- 
fully decomposing the nuclear emission. Although they adopt a 
20% uncertainty of the distance throughout the entire sample, 
the best-estimate distance of NGC 3783 is based on the HST 
photometry as the value quoted above. Nevertheless, Robinson 
et al. (2021) emphasized that the distance of NGC 3783 is quite 
uncertain mainly because the galaxy is rather face-on, leading 
to a large uncertainty of the maximum rotational velocity. The 
strong bar of NGC 3783 may also influence the rotation velocity 
of the galaxy (e.g., Randriamampandry et al. 2015) and affect 
the distance measured from the Tully-Fisher relation. 

Based on the Cosmicflows-3 catalog (CF3, Tully et al. 2016) 
of the distance of ~ 18000 galaxies, Graziani et al. (2019) re- 
construct the smoothed peculiar velocity field, for the first time, 
up to z z 0.05 using a linear density field model. They assume 
a fiducial cosmology (Qm, Qa, Ho) = (0.3,0.7, 75) in the model- 
ing while considering deviations of Ho from the fiducial value, 
so that the reconstructed velocity field does not depend on the 
assumed Hy. They find the cosmic expansion is consistent with 
their fiducial Ho (see also Tully et al. 2016). We use the CF3 
distance-velocity calculator provided by the Extragalactic Dis- 
tance Database (EDD)? to estimate the luminosity distance of 
NGC 3783 based on its vi, corrected for cosmological effects 
(Davis & Scrimgeour 2014; Kourkchi et al. 2020), 


c 


1 1 
vi = Mis |1 5 7 qas - cQ — 0- 3022 |, (9) 


1 
qo = 5 (Om — 2€). (10) 


5 http://edd.ifa.hawaii.edu/. 


Adopting €, = 0.3 and Qa = 0.7, we find qo = —0.55 and 
NGC 3783 has vf, = 2646 km s^! ^ The corresponding distance 
of NGC 3783 is 35.1 Mpc according to the CF3 calculator. It is 
not straightforward to estimate the uncertainty associated with 
this method. We expect that the principal source of uncertainty 
comes from the peculiar velocity of the galaxy due to the non- 
linear effects that cannot be described by the linear model of 
Graziani et al. (2019). They approximate the nonlinear effects in 
the model with a single parameter of nonlinear velocity disper- 
sion onL. Taking the face value of on, « 280km s^! from their 
Bayesian inference, the linear v. of NGC 3783 is in the range 
2366-2926 kms~!. The CF3 calculator yields distances ranging 
from 30.1 to 41.8 Mpc. Therefore, we estimate the uncertainty 
of the distance from the CF3 calculator is at least 6 Mpc.’ 

For completeness, NED provides the estimates of galaxy dis- 
tances based on the multiattractor model by Mould et al. (2000). 
We obtain a Hubble flow velocity 2638 km s^! after correcting 
for the infall due to Virgo cluster, Great Attractor, and Shapley 
supercluster. Taking our fiducial Hy = 70 km s^! Mpc7!, the dis- 
tance of NGC 3783 is 37.9 Mpc, very close to the result of CF3. 
NGC 3783 is found in a group of 9 galaxies according to an up- 
dated nearby galaxy catalog by Kourkchi & Tully (2017). The 
weighted distance using the 2 galaxies with measured distance 
from Cosmicflows-3 is 42.1 + 5.9 Mpc. Although with consid- 
erably large uncertainty, this estimate provides an indication on 
the absolute uncertainty of the indirect approach. 

As displayed in Figure 5, our newly measured distance is 
fully consistent with the results based on other direct and indi- 
rect methods, among which the RMS scatter is ~ 5 Mpc. The 
uncertainty of our measurement is comparable to those of the 
Tully-Fisher relation. 


5.3. Statistical and systematic uncertainties 


The spectral resolution of GRAVITY is moderate and noise of 
adjacent channels is likely correlated. This may lead to an un- 
derestimate of the inferred uncertainty. To gauge the effect of the 
correlated noise, we perform the analysis in two ways, (1) using 
only half of GRAVITY differential phase spectra (every second 
channel) and (2) rebinning across every two channels in the dif- 
ferential phase spectra. The inferred parameter uncertainties are 
typically < 10% higher than those reported in Table B.1. There- 
fore, we conclude that the effect of the correlated phase noise is 
moderate.^ 

The relative uncertainty of our measured distance is about 
33%, which is the combination of the uncertainty of the linear 
radius (Rgır) and the angular radius (Ogjg) of the BLR. We 
find Rgir has a relatively small uncertainty of about 14%, and is 
mainly constrained by the RM data. In contrast, Og; g has about 


^ Itis worth noting that vi, does not depend on Ho (Davis & Scrimgeour 
2014). And the correction with Equation (9) is < 1% as Zis is very small. 

5 Another calculator, NAM, based on nonlinear model of galax- 
ies within 38 Mpc (or 2850 kms!) is provided by EDD. However, 
NGC 3783 is just on the upper boundary of the NAM calculator. It 
yields D; = 37.5 Mpc fully consistent with that of the CF3 calculator, 
although the NAM result is likely much more uncertain as the nonlinear 
model is poorly constrained at the edge of application (Kourkchi et al. 
2020). We therefore prefer the result from the CF3 calculator. 

6 Tn contrast to the 2-o credible interval used in GC20 and GC21, we 
report the 1-c uncertainty level throughout this paper for the simplicity 
of making fair comparisons to the distance measurements with other 
methods. The posterior distributions of most of the key parameters have 
profiles close to Gaussian, so the choice of the credible interval does not 
affect our conclusions. 
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25% uncertainty, mainly due to the relatively large uncertainty of 
the differential phase. The differential phase measurements have 
about 0.1° uncertainty per baseline, so the relative uncertainty of 
the ~ 0.2° phase signal on three baselines is about 30%. There- 


fore, we find 6p, ~ 0S s + ô, where óp, , Óg,, ,, and ôg are the 


fractional uncertainties of the distance, linear radius of the BLR, 
and the differential phase combining 3 baselines. This is overall 
consistent with the conclusion of Songsheng et al. (2021) based 
on tests with mock data. 


A primary concern in terms of systematic uncertainty is that 
the BLR structure of different lines measured by RM and GRAV- 
ITY may be different. A detailed comparison of HG and Bry BLR 
structures for NGC 3783 is out of the scope of this work and 
will be studied in a separate paper where the velocity-resolved 
RM modeling will be used. Alternatively, the broad line profile 
contains the information of the BLR structure. Following Wang 
et al. (2020), we can estimate the relative size difference be- 
tween the H8 and Bry BLRs using (q — 1)/(q + 1), where q = 


(FWHMug JFWHMs;;) . We measure FWHM;, ~ 3680 km s^! 
from our SINFONI spectrum (GC21). As introduced in more de- 
tail in Appendix B.3, we fit globally the H6 complex of 50 RM 
spectra individually, including the narrow and broad components 
of Hf and He i1 44686, [O 111] 444959,5007, and Ferr emission, 
as well as the power-law continuum (e.g., Barth et al. 2015; Hu 
et al. 2015). Accounting for the instrumental broadening (Bentz 
et al. 2021), we find the FWHM of Hf is stable over the RM cam- 
paign and FWHMgg = 4341 + 197 km s7}, Our result is slightly 
smaller than the FWHMmygg of the mean spectrum reported by 
Bentz et al. (2021) mainly due to our different approaches to fit 
the continuum. Therefore, the size difference between the H8 
and Bry BLRs can be about ~ 15%, which is comparable to 
the 1396 RMS of different distance measurements (Section 5.2). 
Calculations with photoionization models will also provide use- 
ful insights (Zhang et al. 2021), although special attention needs 
to be paid to the difficulty of reproducing the observed flux ratios 
of Hydrogen lines in CLOUDY (Netzer 2020). RM observations 
of the same line as that observed with GRAVITY can avoid this 
problem. 


The short time lag of NGC 3783 makes it much more 
straightforward to compare the BLR of different lines from RM 
and GRAVITY than that of 3C 273. The ~100 day time lag 
of 3C 273 requires the campaign to be at least several years. 
The continuum light curve suffers from a further complication 
caused by an overall long-term trend (Li et al. 2020). In fact, 
the newly measured lag is about 2 times smaller than the re- 
sult of early RM campaigns (Kaspi et al. 2000). The dynam- 
ical timescale of the BLR (Peterson 1993) in NGC 3783 is 
tayn ~ Rgrig/VrwuM * 3 years. The BLR structure likely varies 
on a timescale 2 fayn (e.g., Peterson 1993; Peterson et al. 2004; 
Lu et al. 2016). This explains why we find that the Bry spec- 
tra in GRAVITY measurements from 2018 to 2020 show overall 
consistent line profiles, and that these are also consistent with the 
profile from the SINFONI observation in 2019. Although it is not 
possible to entirely exclude it, we do not expect significant size 
variation in the BLR between GRAVITY and RM observations. 
Meanwhile, we emphasize that quasi-simultaneous observation 
of GRAVITY and RM within fayn is necessary to avoid issues 
associated with the variability of the BLR structure. 
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Fig. 5. Comparison of the distances of NGC 3783 measured with the 
methods listed in Table 2. The 1-c uncertainties of direct measurements 
are plotted, while those of indirect measurements are not likely smaller 
than the former. We discuss the uncertainties of the indirect measure- 
ments in the text whenever they can be estimated. 


5.4. Toward a new estimate of the Hubble constant 


With a measurement of D4, one natural step forward is to esti- 
mate the Hubble constant. We derive Ho = vt. / [Dad + aj! = 


65178 km s^! Mpc™! from the joint analysis of NGC 3783. The 
peculiar velocity of NGC 3783 is not included in this estimate. 
Our estimate (Sec. 5.1) suggests the additional systematic uncer- 
tainty introduced by the peculiar velocity is x 10%. At this level, 
the differential phase errors dominate the uncertainty of Hp. 
With upgrades to GRAVITY including adaptive optics with laser 
guide star and wide angle off-axis phase referencing, GRAV- 
ITY+ will allow us to observe at least several hundred suitable 
Type 1 AGN targets up to redshift 2-3, with Pao, Pa, Pay, Ha, 
and Hf redshifted into K band. These lines are much stronger 
than Bry with respect to the continuum. Therefore, we expect 
much less statistical uncertainty (e.g., < 10%) on D4 with future 
observations. Assuming the same sensitivity level as the current 
GRAVITY performance, Songsheng et al. (2021) show that it 
might in principle be possible to obtain a precision measurement 
of Ho using a large sample of differential phase and RM mea- 
surements. In practice, more effort would be necessary to address 
the potential bias due to the properties of different hydrogen lines 
and the assumptions of the BLR model. 


6. Conclusions 


NGC 3783 is the second AGN observed by both GRAVITY in- 
terferometry and RM campaigns. We fitted both data sets simul- 
taneously with a SARM joint analysis. The inferred model pa- 
rameters are fully consistent with the previous study using only 
GRAVITY data, and the uncertainties of key model parameters 
are significantly reduced. In particular, the inferred BH mass is 
2.54710 x 10’ Mo (68% credible interval). For this parameter, 
the uncertainty is a factor of two smaller than both GRAVITY- 
only inference and RM measurements using the integrated Hf 
light curve (dominated by the 0.3 dex calibration uncertainty of 
the virial factor). We also confirm the previous finding of GC21 
that the BLR is highly concentrated with an extended tail of 
clouds out to large radius, which leads to an apparent discrep- 
ancy between the mean radius of the BLR and the measured time 
lag. 
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With the joint analysis, we are able to constrain the angu- 


lar diameter distance of NGC 3783 to be 39007 Mpc. Be- 


cause NGC 3783 is in the nearby Universe, this distance can be 
compared to other independent direct and indirect measurements 
based on the Tully-Fisher relation, galaxy flow models, and its 
galaxy group. The BLR-based geometric distance is fully con- 
sistent with these other results. The dominant uncertainty for the 
distance comes from the differential phase, which has a relative 
uncertainty of about 3096, while the relative uncertainty from 
RM data is about 14%. 

With the ongoing upgrade of GRAVITY, we expect to ob- 
serve many AGNs with improved sensitivity. Our analysis indi- 
cates that the uncertainties of BLR model parameters, such as 
Ogg and Mpgg, will be reduced as the phase uncertainty de- 
creases. Being able to substantially reduce the statistical uncer- 
tainty of Ho looks promising by both improving the precision of 
individual measurements and averaging the measurements over 
many targets. In this context, it will be important to better un- 
derstand the systematic uncertainties associated with the SARM 
method. Further, measuring the same broad line in K band with 
both RM and GRAVITY will be important for future study. 
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Appendix A: The model of the continuum light 
curve 


The observed continuum light curve can be described by (Ry- 
bicki & Press 1992; Li et al. 2018), 


le = s + Eq + ne, (A.1) 


where the vector s is the variable component of the light curve 
with mean zero, Eq represents the mean flux of the light curve, 
and ne is the noise of the observation. E is a vector of unity 
with the same length of the observed light curve and q is the 
mean of the light curve. Assuming that the light curve is a multi- 
variate Gaussian with covariance described by Equation (2) and 
the measurement uncertainties are Gaussian and uncorrelated, 
we can interpolate and extrapolate the light curve with a uni- 
formly and densely sampled continuum light curve model. The 
observed light curve is sampled in a sparse uneven time series, 
tı, meanwhile the model light curve is densely sampled in a uni- 
form time series, t2. The posterior mean and covariance of the 
variable component of the model light curve can be calculated, 


$ = (Cu! C12) (le - Eq), 
C, = C5 - (Ci Ci! Cv, 


(A.2) 
(A.3) 


where Cy, = S (t1, ti)-o tL, Cy. = S (tj, t2), and C2? = S (t5, t2) 
are matrices calculated based on Equation (2), and superscript 
“T” denotes transposition. We denote the observation uncertain- 
ties as o, and the identity matrix as J. The mean and variance of 
the long-term average of the model light curve are, 


ĝ = CE" Ca ly (A.4) 
C, = (E"Cu E". (A.5) 
Therefore, the model continuum light curve is, 

I, = (Lx, + 8) + E( Cx, + d). (A.6) 


where L is the lower triangular matrix of C,, so that C, = LL’, 
x, describes the deviation of the light curve variable component 
from its posterior mean values, and x, describes the deviation of 
the long-term mean from its posterior mean value. We fit x, and 
x4 in the joint analysis as free parameters for the continuum light 
curve model. 


Appendix B: Model fitting results and further tests 


As discussed in Section 2, we use the differential phase data after 
subtracting the continuum phase primarily due to the secondary 
continuum emission component according to the reconstructed 
image (GC21). In order to display the goodness of the fit, we dis- 
play the differential phases of UTA-UT2, UT4-UTI, and UT3- 
UTI, averaged over the three epochs, in Figure B.1 together with 
the best-fit models. The phase signal is dominated by the BLR 
component while the residual continuum phase is moderate. The 
BLR phase signal is not significantly detected in the other three 
short baselines. 

The posterior probability distributions of all of the physical 
parameters are shown in Figure B.2. As discussed in Section 4, 
the posterior distributions of the BLR model of the joint analysis 
show very good consistency with those derived using GRAV- 
ITY data alone (GC21). The uncertainties of several key param- 
eters, such as f, 9, i, and Mpy, are reduced by a factor of ~ 2, 
while those of the other parameters stay the same. The reduc- 
tion of uncertainties of the inferred parameters is also found in 
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Fig. B.1. Averaged differential phase data (points) and the best-fit mod- 
els (lines) of the three baselines that show the strongest signal of the 
BLR component (dashed lines). The best-fit residual continuum phases 
are shown with dotted lines. 


Parameters Model Inference 
Opi g (uas) 11727 
Rgi n (ld) 162575 
Omin (uas) seu 
B post 
6, C) 385 
ic) 30H 
PA (? E of N) 3007 

K —0.03*0.77 
y LI 

é 0574 
Offset (uas) — | (-1.122. 21.622) 
log (Mgu/Ms) 7.40+9-13 
Jenip 0.654913 
P(inflow) 0.79 

6. C) 32:0 
Avgrg (km s^!) 328181 
Da (Mpc) 390905 
log (ra/day) 1.941072 
log (&a/ V day) -0.67 $603 
x 0.672 


Table B.1. Posterior median and 68% credible intervals for the BLR 
model parameters. The central offset of the BLR is in (R. A., Dec.). 
P(inflow) is the probability of inflow where fno, < 0.5 in the posterior 
sample. Avgrr is the difference between the velocity derived from the 
best-fit 44,4, and the systemic velocity based on the redshift. 
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Fig. B.2. Posterior distribution of physical parameters from the Bayesian joint analysis. 


the joint analysis of 3C 273 (Wang et al. 2020). This indicates 
that RM data help to constrain the model in a consistent man- 
ner with GRAVITY data. The B parameter is closer to the upper 
boundary than GRAVITY only inference, reinforcing our previ- 
ous finding that the BLR is heavy-tailed. With GRAVITY data 
only, we found £ peaks at ~ 1 with a long tail to ~ 0, indicat- 
ing that € is not strongly constrained. With the joint analysis, € 
is even less constrained with a nearly flat posterior distribution 
(see Figure B.2). This means that the fitting is not sensitive to the 
midplane obscuration, mainly because most of the BLR clouds 
(65%) are in elliptical orbits rather than radial orbits. The best- 
fit model favors moderate inflow, which is consistent with Bentz 
et al. (2021). 
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Appendix B.1: Constraining the BLR outer radius 


We further test whether the fitting results are significantly differ- 
ent when we constrain the BLR inner and outer radii based on 
previous RM of high-ionization lines and NIR continuum. Fol- 
lowing the discussion of GC21, we cast additional prior informa- 
tion that Rmin = 4 ld and Rout = 80 Id. The key results are shown 
in Table B.2, where D4 is reduced by about 30%. The angular 
size of the BLR is almost the same as our primary model without 
additional constraints, as it is directly constrained by the GRAV- 
ITY phase data. At the same time, Rp; g is reduced by about 30% 
due to the hard boundary Roy = 80 ld of the cloud distribution. 
Given the observed angular size of the BLR, a smaller D4 is re- 
quired in order to incorporate a smaller BLR as required by the 
prior. In addition, the inferred 6 is around 1, meaning that the 
cloud distribution is now exponential rather than heavy-tailed. 


GRAVITY Collaboration: A geometric distance to the supermassive black Hole of NGC 3783 


This is because Rout limits the extent of the cloud distribution. 
The peak of the model differential phase is correspondingly re- 
duced, although the model fits the data almost as well as our 
primary model given the data uncertainty. All of the other in- 
ferred model parameters are consistent with our primary model 
results. Nevertheless, the radius-constrained model is not pre- 
ferred over our primary model because the inferred D4 would be 
the smallest of the measurements discussed in Section 5. Quan- 
titatively, it is about 1.5-c lower than the other measurements, 
which cautions against the additional constraints to the model. 
While Rout = 80 ld is a reasonable choice based on NIR contin- 
uum RM analyses (Glass 1992; Lira et al. 2011), the contribution 
of NIR emission from the accretion disk may bias the time lag to 
a small value (Lira et al. 2011). Nevertheless, the joint analysis 
may also be biased by the assumptions of the BLR model. More 
comprehensive comparisons will be performed in future paper 
that considers different model assumptions, e.g., different cloud 
radial distributions and BLR size constraints. 


Appendix B.2: The nonlinear response of RM 


The BLR reverberates the Hydrogen ionizing (UV) photon vari- 
ability from the accretion disk. The optical (i.e., V band) light 
curve is usually used as an approximation to the UV light curve. 
However, their correlation is not necessarily linear. Following 
(Li et al. 2018, see also Li et al. 2013), we revise Equation (3) 
into 


lug(t) = A J P(r) 1*9 (t — 7) dr, (B.1) 


where 6, captures the nonlinearity of the response. The model 
inference is conducted in the same way for the other parameters 
as discussed in Section 3, with the uniform prior of 6, between 
-1 and 3 (Li et al. 2018). We find the inferred model parame- 
ters are entirely consistent with our primary fitting results (Ta- 
ble B.1). Some of the key parameters are listed in Table B.2. 
Meanwhile, the uncertainty of Rgrg increases by about 70%, 
resulting in an increase to the uncertainty of D4 by the same 
amount. This mainly originates from the degeneracy between 6, 
and the scaling factor A (Equation (B.1)). Optimizing the sam- 
pling method will reduce the uncertainties of some parameters 
due to the degeneracy. However, this is beyond the scope of the 
current paper. On the other hand, the fitting of the light curve is 
not substantially improved. Therefore, we opt not to include the 
nonlinear effect of the RM response in the modeling. 


Appendix B.3: Measuring the light curve using the 
decomposition method 


Here we test whether the joint analysis is sensitive to the method 
measurement the light curve. Bentz et al. (2021) measured the 
Hf flux by fitting a local linear continuum underneath the Hg 
line. This method has been widely used and has been shown to 
be effective in producing accurate light curves for strong emis- 
sion lines, such as Hf (e.g., Peterson et al. 1998; Kaspi et al. 
2000; Grier et al. 2012). Recent studies have developed methods 
to decompose the emission lines by simultaneously fitting vari- 
ous physical components to the spectra over a wide wavelength 
range (Barth et al. 2015; Hu et al. 2015; Shen et al. 2016). We 
therefore conducted a joint fit using the Hf light curve measured 
with the decomposition method in order to test whether the joint 
analysis is sensitive to the measurement of the light curve. Fol- 
lowing the approach of Barth et al. (2015), our model consists of 


the AGN power-law continuum, [O i11] doublet, broad and nar- 
row Hf and Het lines, and the Ferr blended lines. We find the 
stellar continuum, reddening, and Her lines at 4471, 4922, and 
5016 À are not necessary to achieve a good fit so, to improve 
the stability of the decomposition of different spectra, we do not 
include these components. The broad and narrow emission lines 
are all fitted with a fourth-order Gauss-Hermite function (van 
der Marel & Franx 1993), whose centroid, dispersion, h3, ha, 
and amplitude are the parameters in the fitting. The profiles of 
the other narrow lines (the dispersion, h3, and h4) are tied to 
that of [O 111] 45007. The ratio of [O 111] 444959,5007 is fixed 
to 1/3. The rest-frame wavelength separation of the two lines is 
fixed to 47.9 A. The Fer emission lines are modeled with the 
template from Boroson & Green (1992). We adjust the redshift, 
dispersion, and amplitude of the template in the fitting. We fit 
the spectra at rest-frame 4500—5500 A, masking the small wave- 
length range of [N 1] lines at 5199 and 5201 A. We subtract all of 
the best-fit components except the narrow and broad Hf line, and 
the H£ flux is integrated over rest-frame 4800—4940 A, in order 
to be consistent with Bentz et al. (2021). To estimate the uncer- 
tainty of the light curve, we perturb the data according to the 
measurement uncertainty of the spectra and redo the measure- 
ment. The final uncertainty is the quadrature sum of the statisti- 
cal uncertainty from the bootstrapping and 1% of the integrated 
flux to make sure that the resulting uncertainty is not too small. 

The light curve measured by the decomposition method has 
a very similar shape to that from Bentz et al. (2021) as shown 
in Figure B.3. The former is about 5% higher than the latter, 
probably because the best-fit continuum from the decomposi- 
tion method is slightly lower than that from the local continuum 
underneath the Hf line. Using the method introduced in Sec- 
tion 3, we are able to infer model parameters consistent with 
those reported in Table B.1 within 1-c-. However, the peak-to- 
valley variation amplitude of the new light curve is reduced by 
about 2096. Correspondingly, the new inference prefers about 
30-50% larger BLR radius, inclination angle, and D4 in order to 
reduce the variation amplitude of H£ light curve. The model in- 
ference is slightly biased in this way. To compensate for this, we 
need to take the nonlinearity of the H8 response (Equation B.1) 
into account. This yields inferred parameters almost exactly the 
same as those reported in Table B.1 with the additional param- 
eter óc ~ —0.2. The key parameters are listed in Table B.2. The 
responsivity of H£ is physically expected to be x 1 (i.e., 6. < 0) 
based on photoionization models (e.g., Gaskell & Sparke 1986; 
Korista & Goad 2004) as well as observations of e.g., NGC 5548 
(Goad & Korista 2014). We emphasize that the observational ef- 
fects, such as using the V-band continuum as a proxy for the 
ionization continuum and spectral decomposition over a wide 
wavelength range, may also introduce a nonlinear effect. The 
change of the light curve variation amplitude does not influence 
much the time lag of the RM. However, our test shows that the 
joint analysis may be biased if the nonlinear effect is not properly 
taken into account. 
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Model Inference with 
Parameters 
(1) Boundary (2) Nonlinearity (3) New light curve 
Oni (Has) 71413 7047 7471 
Rgi n (ld) 12e 18:935 16.4*63 
D4 (Mpc) 30.0775 43.6277 36. 8 
log (Mpu/Mo) Tole 7.407045 7.501923 


Table B.2. Same as Table 1, but showing the fitting results: (1) when Rmin = 4 ld and Rmax = 80 ld are required (Appendix B.1); (2) when the 
nonlinear response is considered (Appendix B.2); (3) when we use the light curve measured by the decomposition method (Appendix B.3). 
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Fig. B.3. Comparison of the H£ light curves measured by Bentz et al. 
(2021, in blue) and the decomposition method (in red, see the text). 
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